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Abstract. We present an algorithm which takes as input a closed semi- 
algebraic set, S C R fc , defined by 

Pi < 0, . . . , P e < 0, Pi G R[X 1} X k ], deg(Pi) < 2, 



and computes the Euler-Poincare characteristic of S". The complexity of 
the algorithm is . 
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\q ! 1. Introduction 

O 

Let R be a real closed field and let 5* C R be a basic semi- algebraic set defined 
by P 1 <0,...,P e < 0, with P t G R[X 1 ,...,X k ),deg{Pi) < 2, 1 < % < i. 
It is known [21 E] that the sum of the Betti numbers of S (and hence also its 
Euler-Poincare characteristic) is bounded by k 0( - e \ Notice that for fixed £, these 
bounds are polynomial in k. One can also check whether S is non-empty, as well 
as compute a finite set of sample points meeting every connected component 
of 5* in time However, no algorithm with similar complexity is 

known for computing any of the individual Betti numbers of 5* (for instance, the 
number of connected components). The best known algorithm for computing 
all the Betti numbers of S has complexity k 2 

Here, and elsewhere in this paper the Betti number, bi(S), is the dimension 
of the simplicial homology group, Hi(S,Q), in case S is closed and bounded. 
If 5 is a closed but not necessarily bounded semi-algebraic set, bi(S) is the 
dimension of H^S D 5^(0, r),Q), for sufficiently large r > (here and in the 
rest of the paper, Bk(0,r) denotes the open ball of radius r in R h centered at 
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the origin, and X denotes the closure of a semi-algebraic set X). It is easy to 
see that bi(S) is well-defined and we denote by 

x(s) = £(-1)^0?) 

8=0 

the Euler-Poincare characteristic of S. 

In this paper we describe an algorithm for computing the Euler-Poincare 
characteristic of S, whose complexity is k°^\ Our algorithm relies on an 
efficient algorithm for computing the Euler-Poincare characteristic of the real- 
izations of all realizable sign conditions of a family of polynomials described in 
and uses different techniques than those used in |2J [TUJ . 

The main result of this paper is the following. 
Main Result: We present an algorithm Q Algorithm 4.2 in ISection 4|) which 



given a set of I polynomials, V = {Pi, . . . , Pg} C R[Xi, . . . , Xk], with deg(p) < 
2, 1 < i < £, computes the Euler-Poincare characteristic, x(S), where S is the 
set defined by Pi < 0, . . . , Pi < 0. The complexity of the algorithm is . 
If the coefficients of the polynomials in V are integers of bitsizes bounded by 
r, then the bitsizes of the integers appearing in the intermediate computations 
and the output are bounded by rk°( e \ 

The rest of the paper is organized as follows. In ISection 2\ we describe some 
mathematical and algorithmic results we will need for our algorithm. We also 
include a brief introduction to spectral sequences since they play a motivating 
role in the design of the main algorithm described in this paper. In ISection 31 
we describe an algorithm for computing the Euler-Poincare characteristic of 
a set defined by homogeneous quadratic inequalities. Finally, in ISection 41 we 
describe our algorithm for the general (inhomogeneous) case. 

For referring to well known results in real algebraic geometry we sometime 
use reference |7j as a useful source. 



2. Preliminaries 

In this section we describe some mathematical and algorithmic results we will 
require in the rest of the paper. 

2.1. Definition of the Euler-Poincare Characteristic. For the purposes 
of our algorithm, it will be useful to define Euler-Poincare characteristic for 
locally closed semi-algebraic sets. We do this in terms of the Borel-Moore 
homology groups of such sets (defined below). This definition agrees with the 
definition of Euler-Poincare characteristic stated earlier for closed and bounded 



Computing Euler-Poincare Characteristics 



3 



semi- algebraic sets. They may be distinct for semi-algebraic sets which are 
closed but not bounded. 

The simplicial homology groups of a pair of closed and bounded semi- 
algebraic sets T C S C R fc are defined as follows. Such a pair of closed, 
bounded semi-algebraic sets can be triangulated [7j using a pair of simplicial 
complexes (K, A), where A is a sub-complex of K. The p-th simplicial homol- 
ogy group of the pair (S, T), H p (S, T), is H P (K, A). The dimension of H p (S, T) 
as a Q-vector space is called the p-th Betti number of the pair (S, T) and 
denoted b p (S,T). The Euler-Poincare characteristic of the pair (S,T) is 

x(s,t) = "£(-i)%(s, T). 

i 

The p-th Borel- Moore homology group of S C R fc , denoted H^ M {S), is 
defined in terms of the homology groups of a pair of closed and bounded semi- 
algebraic sets as follows. For r > 0, let S r = SnB k (0, r). Note that, for a locally 
closed semi- algebraic set S, both S r and S r \ S r are closed and bounded and 
hence H p (S r , S r \ S r ) is well defined. Moreover, it is a consequence of Hardt's 
triviality theorem [TT] that the homology group H p (S r , S r \ S r ) is invariant 
for all sufficiently large r > 0. We define, H^iS) = H p (S~ r ,S~ r \ S r ), for 
r > sufficiently large, and it follows from the remark above that it is well 
defined. The Borel-Moore homology groups are invariant under semi-algebraic 
homeomorphisms (see It also follows clearly from the definition that for 
a closed and bounded semi- algebraic set, the Borel-Moore homology groups 
coincide with the simplicial homology groups. 

For a locally closed semi-algebraic set S, we define the Borel-Moore Euler- 
Poincare characteristic by, 

x BM (s) = J2bF M (s), 

where bf M (S) denotes the dimension of H^ M (S, Q). Note that if S is closed 
and bounded, then x BM (S) = x(S). 

The Borel-Moore Euler-Poincare characteristic has the following additive 
property. 

Proposition 2.1. Let X, Xi and X 2 be locally closed semi-algebraic sets such 
that 

iiux 2 = x, x 1 nx 2 = 0. 

Then 

x BM (x) = x BM (x 1 ) + x BM (x 2 ). 
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Proof. This is classical (see for example, Proposition 2.6 in [U| for a proof). 

□ 

2.2. Sign Conditions and their realizations. A sign condition is an ele- 
ment of {0, 1, —1}. We define 

{0 if and only if x = 
1 if and only if x > 
— 1 if and only if x < 

Let Z C R fc be a locally closed semi-algebraic set and let V = {Pi, . . . , P s } 
be a finite subset of R[Xi, . . . ,Xk]. A sign condition on V is an element of 
{0,1,-1}*. 

The realization of the sign condition a on Z is 

ft(<7, Z) = {x e Z | f\ sign(P(x)) = <r(P)}, 

and its Euler-Poincare characteristic is denoted y^ M {o, Z). 

We denote by Sign(P, Z) the list of a e {0,1, -1}^ such that K(a, Z) is 
non-empty. We denote by x™^, Z) the list of Euler-Poincare characteristics 
X BM (a, Z) = X BM (n^ Z)) for a G Sign(P, Z). 

We will use the following algorithm for computing the list x SM (^ > ) Z) de- 
scribed in 0. We describe here the input, output and complexity of the algo- 
rithm. 

Algorithm 2.2. Euler-Poincare Characteristic of Sign Conditions. 

Input: an algebraic set Z = Z(Q, R fc ) C R fc and a finite list V — Pi, . . . , P s of 

polynomials in R[Ai, . . . , X^]. 
Output: the list X BM (V,Z). 

Complexity: Let k' be the dimension of Z, d a bound on the degree of Q and 
the elements of V and s = #(V)). The number of arithmetic operations is 

s k ' +1 0(d) k + s k '((k'\og 2 (s) + k\og 2 {d))d)° [k) ■ 

The algorithm also involves the inversion matrices of size s k '0(d) k with integer 
coefficients. 

If D = Z, and the bitsizes of the coefficients of the polynomials are bounded 
by r, then the bitsizes of the integers appearing in the intermediate computa- 
tions and the output are bounded by r((/c'log 2 (s) + k\og 2 (d))d)°^ k \ 
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2.3. Infinitesimals. In our algorithms we will use infinitesimals in order to 
ensure that the set we are dealing with is bounded. To ensure this we will 
extend the ground field R to R{e), the real closed field of algebraic Puiseux 
series in e with coefficients in R. The sign of a Puiseux series in R(e) agrees 
with the sign of the coefficient of the lowest degree term in e. This induces 
a unique order on R(e) which makes e infinitesimal: e is positive and smaller 
than any positive element of R. 

If R' is a real closed field containing R, then given a semi-algebraic set S 
in R fc , we denote the extension of S to R' by Ext(S', R'). Ext(S f , R') is the 

fk 

semi- algebraic subset of R defined by the same quantifier free formula that 
defines S. The set Ext (5", R') is well defined (i.e. it only depends on the set 
S and not on the quantifier free formula chosen to describe it). This is an 
easy consequence of the Tarski-Seidenberg transfer principle (see for example 
Section 2.4.1 in 0). 

2.4. Spectral Sequences. For the benefit of the readers we include a brief 
introduction to the theory of spectral sequences pointing to [UJ H2] for more 
details. 

A spectral sequence is a sequence of bigraded complexes (E r , d r : E™ — ► 
£P+r,q-r+i^ suc ] 1 that the complex E T+ \ is obtained from E r by taking its ho- 
mology with respect to d r (that is E r+ i = H dr (E r )). 

There are two spectral sequences, 'E^ q , "El' q , (corresponding to taking row- 
wise or column-wise nitrations respectively) associated with a double complex 
C"'*, which will be important for us. Both of these converge to #*(Tot*(C"'*)). 
This means that the homomorphisms d r are eventually zero, and hence the 
spectral sequences stabilize, and 

'E™ = "E™ - iP(Tof(C"'-)), 

p+q=i p+q=i 

for each % > 0. 

The first terms of these are 'E x = H S (C*>'),'E 2 = H d H 5 (C'''), and "E 1 = 
H d (C"),"E 2 = H s H d (C'>-). 

In particular, assuming that the complex C*'' is bounded in both directions, 
we have that, 

Proposition 2.3. 

E i>0 (-l)Mim(^(Tof(C-'-))) = £ M>0 (-irMim('i^) 

= E P ^o(-l) P+,? dimr^). 
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2.5. Leray Spectral Sequence of a map. The Leray spectral sequence of 
a proper map, 

is a classical tool in algebraic topology which relates the cohomology of the 
space A, with those of the space B and of the fibers of the map /. Its most 
common use is in the theory of sheaf cohomology jH]- We will need it in a 
very special situation where the sets A and B are compact semi-algebraic sets 
and / a continuous semi-algebraic map. In this special situation it is possible 
to define the Leray spectral sequence in terms of triangulations, which we do 
below. 

Consider a semi- algebraic continuous map, / : A — > B, where A and B 
are compact semi-algebraic sets. Moreover, let h : A — > B be a semi-algebraic 
triangulation of B, and let Tl(A) denote a cell-complex, such that A is the union 
of the cells in TC(A), and such that for any simplex a G A, A a = / _1 (/i(er)) is 
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a subcomplex of TC(A) ( where X denotes the topological closure of X). 

Then, the Leray spectral sequence of / is isomorphic to the spectral sequence 
(corresponding to the column-wise filtration) associated to the double complex 
C°'* defined as follows: 

C**= C%A a ), 

a a p-simplex in A 

where C q {A a ) denotes the vector space of g-co-chains of the complex A a . The 
horizontal and the vertical differentials are the obvious ones (see 0). The 
spectral sequence associated to the double complex defined above converges to 
the co-homology of A. 

3. The basic homogeneous case 

Let V = {Pi, . . . , Pi} C R[X , Xi, . . . , X k ] be a set of homogeneous quadratic 
polynomials, and let S be the basic closed semi- algebraic set defined on the 
unit sphere S k C R fc+1 by the inequalities, 

Pi < 0, . . . , Pi < 0. 

We denote by Si the subset of S k defined by Pi < 0. Then, S = C\ e i=1 Si. 
For J C { we denote by S J = Uj e jSj. The following equality is a 
consequence of the exactness of the Mayer- Vietoris sequence. 

Lemma 3.1. 

X(S)= J2 (-!) #(J)+1 X(^)- 

JC{1,...,£} 

Proof. In case £ = 2, this is a direct consequence of the exactness of Mayer- 
Vietoris sequence (see for example, Corollary 6.28). The general case follows 
from an easy induction. □ 

Thus, in order to compute x{S) ^ suffices to compute x(S J ) for each J C 
{1,. ..,£}. 

3.1. Topology of sets defined by quadratic constraints. We first recall 
some facts about topology of sets defined by quadratic inequalities [T]. Let 
Pi, . . . , P s be homogeneous quadratic polynomials in R[X , . . . , Xk]- 

We denote by P = (Pi,..., P.) : R fe+1 -> R s , the map defined by the 
polynomials P\, . . . ,P S . Let 

A = [J {x G S k I Pi(x) < 0}. 

1<1<S 
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Let 



Q= {cue R s \ \u\ = 1, Ui < 0, 1 < % < s}. 



For oo G f2 we denote by ujP the quadratic form denned by 



Let B C £1 x S k be the set defined by, 



B = {(u, x) | u G Q, x G S k and ljP(x) > 0}. 



We denote by 0i '■ B —* Q and 02 : B — > £> the two projection maps. 



The following was proved by Agrachev jj. With the notation developed 
above, 

Proposition 3.2. Tie map 02 gives a homotopy equivalence between B and 



Proof. We first prove that <p2{B) — A. If x G A, then there exists some 
i, 1 < i < s, such that Pj(x) < 0. Then for u; = (— 5u, . . . , — 6 si ) (where 5^ = 1 
if i = j, and otherwise), we see that (co,x) G B. Conversely, if x G 4>2(B), 
then there exists c<j = (u>i, . . . , u> s ) G such that, J2t=i UiPi(x) > 0. Since, 
^<0,l<i<s, and not all uj{ = 0, this implies that Pj(x) < for some 
i,l < i < s. This shows that x E A. 
For x G 02 {B), the fibre 
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MB) = A. 



2 (x) = {(u),x) \u eQ such that ujP(x) > 0} 



can be identified with a non-empty subset of defined by a single linear in- 
equality. From convexity considerations, all such fibres can clearly be retracted 
to their center of mass continuously, proving the first half of the proposition. □ 
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For any quadratic form Q, we will denote by index(Q), the number of 
negative eigenvalues of the symmetric matrix of the corresponding bilinear 
form, that is of the matrix M such that, Q(x) = (Mx,x) for all x G R k+1 . We 
will also denote by Aj(Q),0 < i < k, the eigenvalues of Q, in non- decreasing 
order, i.e. 

Ao(Q)<Ai(Q)<--.<A fc (Q). 

Given a quadratic map P = (Pi, . . . , P s ) : R fc+1 — > R s , and < j < k, we 
denote by 

Qj = {ueQ I Xj(ujP) > 0}. 

For notational convenience, will denote the empty set and Qk+i the whole 
space Q. 

It is clear that the Q^s induce a filtration of the space Q, i.e., Qq C Qi C 
■ ■ ■ C fifc+i- 

Agrachev PQ showed that the Leray spectral sequence of the map 0i (con- 
verging to the cohomology H*(B) = H*(A)), has as its Ei terms, 

(3.3) £| 5 = /r p (n fc _g,n fc _ ff _ 1 ). 

This follows from the fact that the fibre of the map 0i over a point ui G 
Qj \ ilj-i has the homotopy type of a sphere of dimension k — j. To see this 
notice that for uj G Qj X (u>P), . . . , A 3 -_i(o;P) < 0. Moreover, letting 

Yq(uiP), . . . , Yk(uP) be an orthonormal basis consisting of the eigenvectors of 
tuP, we have that ^^(cu) is the subset of S k defined by, 

k 

k 

Since, Aj(a;P) < 0,0 < i < j, it follows that for uj G f2j \ 0f 1 (u;) 
is homotopy equivalent to the (k — j) -dimensional sphere defined by setting 
Y (ujP) = ■■■ = Yj^iujP) = on the sphere defined by £*L Y^uP) 2 = 1. 

The following proposition relates the Euler-Poincare characteristic of the 
set A with those of Qj \ Qj-i, < j < k + 1. 

Proposition 3.4. 

fc+i 

X (A) = X BM {A) = $> BM (% \ %_x)(l + (-l)( fe ^). 

j=0 
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Proof. Notice that the sets flj \ are locally closed, and the fibre over 
a point uo G Qj \ Qj-i is compact and has the homotopy type of a (k — j)- 
dimensional sphere. Now consider a sufficiently fine triangulation of Q, which 
respects the filtration fi C ••• C &>k+i, and such that over each simplex a of 
the triangulation lying in Qj (j>x (cr) is homotopy equivalent toax S k ~i . 

The Euler-Poincare characteristic of a (k— j)-dimensional sphere, S k ~\ is equal 
to 1 + (— l)^ - ^. The proposition now follows from the additivity property of 
the Borel-Moore Euler-Poincare characteristic and Proposition 3.2| 



□ 



Since [Proposition 3.4| is central to the algorithm presented in this paper, we 
include a different proof below which uses the spectral sequence First 
note that by [Proposition 2~ 



^(n i \n i _ 1 ) = x (n i )-x(n i _i). 



It follows from the convergence of the spectral sequence in |(3.3)| and Proposition 2.3 
that, 

x(A) = E P+g=i (-i) i dim(£;r) 

= Ep+^-iydimif^^lWi) 

= Eo< 9 < fc+ i Eo< P <k(-l) P+q dim(H^ q , iwO) 

= Eo< ? < fc+ i(- 1 J 9 Eo<p< fc (-l) P dim(^(O fe _ g ,fi fe _ ? _ 1 )) 

Now, from the exact sequence of the pair (Qk-q, &>k-q~i), namely, 



we get that, 



^(-l) l (dim(if l (fi fc _,_ 1 )) - dim^Cfifc.,)) + dim(# l (lV 9 , fi*- g -i)) = 0, 

j>0 

which yields 

0<p<fc 
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Thus, the previous sum 

= E < q < k+ x(-mx(n k - q ) - x(n*-*-i)) 

= Eo< J<fc+ i(-i) fc+1 ^(x(^--i) - x(%- 2 )) 

= Eo< J<k+ i(-^ k+1 - J x(^-i) - E < j<k+ x(-^ k+1 - j x(^-2) 

= Eo^-^xm - Eo<i< fc -i(-i) A:+1 - i x(^) 

= EoW-i)*-^) + Eo^^fc-iC-i)*-^) 

= x (fi fc ) - 2 X (fi fe _i) + 2 X (fi fc _ 2 ) + • ■ ■ + (-i) fc 2 X (n ) 

= Eo< j<fc+ i(x(^) - + (-i)*-') 

= Eoi< fc+ iX BM (^\^-i)(i + (-i) fc ^)- 

□ 

Let Z = (Zi, . . . , Z s ) be variables and let M(Z) be the symmetric matrix 
corresponding to the quadratic form Z ■ P — Z\Px + • — h Z S P S . The entries of 
M(Z) depend linearly on Z. Let 

F(Z, T) = det(M(Z) + T • I k+1 ) = T k+1 + C k T k + ■ ■ ■ + C , 

where each Cj G R[Zi, . . . , Z s ] is a polynomial of degree at most k + 1. 

It follows from the well known Descarte's rule of signs (see for example, 
Remark 2.42 in [7]) that for any z G f2, index(zP) is equal to the number of 
sign variations in the sequence C (z), . . . ,Ck(z),+l. Thus, the signs of the 
polynomials C , . . . , determine the index of zP. For a G {0, +1,— 1} C a 
sign condition on the family C = {C , . . . , C^}, let n(a) denote the number of 
sign variations in the sequence, cr(C ), . . . , cr(Cfc), +1. Let Sign(C, Q) be the set 
of sign conditions realized by the family C on Q. The following proposition is 
an immediate consequence of [Proposition 3.4 and the additivity of the Euler- 
Poincare characteristic. 

Proposition 3.5. 

X(A) = X BM {A) = Yl X BM (K(a, fi)) • (1 + (-1)<*^<*»). 
o-eSign(c,n) 

Before proceeding further we discuss a small example. 

Example 3.6. Let I = 2, k = 2 and P = [P u P 2 ) : R 3 -> R 2 be the quadratic 
map with, 

p v2 I -ia2 x^2 

H — A + A l — A 2) 
p _ y2 y2 Y 2 
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In this example, 

Q — {(ui,U 2 ) | uj 2 + = 1,uj 1 ,uj 2 < 0} 

consists of the arc of the unit circle in the third quadrant. 
Also, 

A = {x <E S 2 \ Pi(x) < V P 2 (x) < 0} 

= {(x , x u x 2 ) e R 3 I x 2 + x\ + x\ = 1, x\ + x\ > 1/2}. 

The set A in this example clearly has the homotopy type of a circle, and hence, 

(3.7) X (A) = 0. 

Now, for uj = (cui,U! 2 ) G f2, 

UoP = UOiPi + LU 2 P 2 

= (ui + u 2 )Xl + (ui - cu 2 )Xf - {ui + uj 2 )X\. 

Following notations introduced above, 

F(Z 1 , Z 2 , T) = (Z 1 + Z 2 + T){Z X -Z 2 + T)(-Z 1 -Z 2 + T) 
= T 3 + (Zi - Z 2 )T 2 - (Zi + Z 2 ) 2 T + (Zi + Z 2 ) 2 (Z 2 - Zi). 

Thus, in this example, 

C (Z 1 ,Z 2 ) = (Z 1 + Z 2 ) 2 (Z 2 -Z 1 ), 
C 1 (Z 1 ,Z 2 ) = -(Z 1 + Z 2 ) 2 , 
C 2 (Zi, Z 2 ) — Z\ — z 2 . 

There are three realizable sign conditions on C = {C , Ci, C 2 , +1} on f2. 
They are, 

o-i = (-,-,+,+), 

<7 2 = (0,-,0,+), 
^3 = ( + ,-,-,+)• 

We have, 

71(01) = 1, 
71(0-2) = 1, 
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n(a 3 ) = 2. 

The realizations TZ(ai,Q) and TZ(a 3 ,Q) are each homeomorphic to [0,1) 
while 7£(cr 2 ,f2) is a point. Thus, 

x BM (a 1 ,n) = x BM (a 3 ,n)=o, 

while 

Finally, we have 

EUx BM (^m-(-^ 2 - n{aj) ) = o(i + (-i) i ) + i(i+(-i) 1 ) 

+ o(i + (-i) 2 ) 
= o, 

which agrees with pT)| 

We are now in a position to describe an algorithm for computing the 
Euler-Poincare characteristic of a union of sets, each defined by a homoge- 
neous quadratic inequality. In the algorithm we will use the following notation. 
Given two finite families of polynomials, V C V, and o G {0, +1, — u' G 
{0, +1, -if, we say that a -< o' iff for all P G P, a(P) = a'{P). 

Algorithm 3.8. Euler-Poincare characteristic of a union. 
Input: A set of quadratic forms {Pi, . . . , P s } C R[X , . . . , Xk]. 
Output: x{A), where A is the set defined on the unit sphere S k C R fc+1 by the 
formula 

Pi < V ■ • ■ V P e < 0. 

1. Let P = (Pi, . . . , P s ). Let Z = (Z u . . . , Z s ) be variables and let M(Z) be 
the symmetric matrix corresponding to the quadratic form Z ■ P. Compute 
the polynomials Cj G R[Zi, . . . ,Z S ] by computing the following determi- 
nant. 

det(M(Z) + T ■ I k ) = T k+1 + C k T k + • • • + C . 

2. Compute x BM {C,Q) as follows. Call [Algorithm 2.2| with input C = C U 
{Zi, . . . , Z s } and Q — Z\ H — • + Z^ — 1. The output is the list 

X BM (C',Z(Q,R k )). 

For each a G {0, +1, -1} C , such that there exists a' G Sign(C, Z(Q, R k )) 
with a -< a' and u'(Zj) G {0, —1} for 1 < j < s, compute 

X BM (cr,Q)= Yl X BM (a',Z(Q,R k )). 

ct' ,a~<a' ,a'(Zj)£{0,-l},l<j<s 
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3. Output 

X (A) = X BM {n{a^)) ■ (1 + (-!)(*-«(-))). 

o-eSign(c,n) 



Proof of Correctness: The correctness of the algorithm is a consequence 
of Proposition 3.5 and the correctness of [Algorithm 2.2 □ 
Complexity Analysis: The complexity of the algorithm is using the 
complexity of Algorithm 2.2 □ 



We are now in a position to describe the algorithm for computing the Euler- 
Poincare characteristic in the basic, homogeneous case. 

Algorithm 3.9. The basic homogeneous case. 

Input: A set of quadratic forms {Pi, . . . , Pi} C R[X , . . . , X^]. 
Output: x(S)i where S is the set defined on the unit sphere S k C R fc+1 by the 
inequalities, 

Pi < 0, . . . , P e < 0. 

1. For each subset J C {1, . . . ,£} do the following. 

2. Compute using Algorithm 3.8\ x{S J ). 

3. Output 

x(S)= Yl (- 1 )* {J)+1 x(s J ). 

Jc{l,...,£} 



Proof of Correctness: The correctness of the algorithm is a consequence 
of Lemma ILemma 3.11 and the correctness of Algorithm 3.8| □ 
Complexity Analysis: There are 2 e calls to [Algorithm 3.8| Using the com- 
plexity analysis of Algorithm 3.8 , the complexity of the algorithm is bounded 
byk°®. □ 



4. The General Case 

Let V = {Pi, . . . , P e } C R[Xi, . . . , X h \ with deg(P i ) < 2, 1 < % < 
S C R fc be the basic semi-algebraic set defined by Pi < 0, . . . , Pi 
< £ be an infinitesimal, and let 



£, and let 
< 0. Let 



P 



3=1 
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Let S' C R(e) k be the set defined by P x < 0, . . . , P i+1 < 0. 

Denoting by P/ 1 the homogenization of Pj, and S C S h the set defined by, 
P[ l < 0, . . . , Pg < 0, < 0, on the unit sphere in R(e) k+1 we have, 

Proposition 4.1. For < z < k, x{S) = x(S') = ~x{S h ). 

Proof. Using the well known conic structure at infinity of semi-algebraic 
sets (see for example Proposition 5.50, [Zj) we have that for all sufficiently 
large r > 0, S nP>k(0,r) is a semi-algebraic deformation retract of S. Since e is 
an infinitesimal, it follows that 5" = Ext(S', R(e)) nPfc(0, -) is a semi-algebraic 
deformation retract of Ext (S, ~R(s)). This implies that y(5') = y(5"). 

To prove the second equality, first observe that S' is bounded, and S h is the 
projection from the origin of the set 1 x S' C 1 x K(e) k onto the unit sphere in 
R fc+1 . Since, S' is bounded, the projection does not intersect the equator and 
consists of two disjoint copies in the upper and lower hemispheres, and each 
copy is homeomorphic to S'. □ 

Algorithm 4.2. The general case. 

Input: A family of polynomials V = {Pi, . . . , Pi} C R[Xi . . . , X k ], with deg(Pj) < 
2. 

Output: x{S)i where 5* is the set defined by 

S = p| { x g R fc I P{x) < 0}. 
Pev 



1. Repla ce the family V by the family, V h = {P/ 1 , . . . , Pf, P e h +1 }. 

2. Using [Algorithm 3. 9| compute x(S h ). 

3. Output X (S) = lx(S h ). ' 



Proof of Correctness: The correctness of [Algorithm 4.2| is a consequence 



of Proposition 4.1 and the correctness of [Algorithm 3.9[ □ 



Complexity Analysis: The complexity of the algorithm is clearly from 
complexity analysis of Algorithm 3.9[ □ 



Remark 4.3. In this paper we have described an algorithm for computing the 
Euler-Poincare characteristic of a basic closed semi-algebraic set defined by a 
constant number of quadratic inequalities P < 0, P G V. It is straightforward 
to extend the algorithm to the case of semi-algebraic sets defined by Boolean 
formulas without negations, whose atoms are of the form P > or P < for 
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P EV, using the technique used in jl] for reducing the problem of computing 
Euler-Poincare characteristic of such sets, to the basic closed case. This reduc- 
tion works perfectly well even in the quadratic situation and does not worsen 
the complexity. 

Acknowledgements 

The author was supported in part by an NSF Career Award 0133597 and a 
Sloan Foundation Fellowship. 

References 

[1] A. A. Agrachev, Topology of quadratic maps and Hessians of smooth 
maps, Algebra, Topology, Geometry, Vol 26 (Russian), 85-124, 162, Itogi 
Nauki i Tekhniki, Akad. Nauk SSSR, Vsesoyuz. Inst. Nauchn.i Tekhn. In- 
form., Moscow, 1988. Translated in J. Soviet Mathematics. 49 (1990), no. 
3, 990-1013. 

[2] A.I. Barvinok, Feasability Testing for Systems of Real Quadratic Equa- 
tions, Discrete and Computational Geometry, 10:1-13 (1993). 

[3] S. Basu, On different bounds on different Betti numbers, Discrete and 
Computational Geometry, Vol 30, No. 1 (2003). 

[4] S. Basu, On Bounding the Betti Numbers and Computing the Euler Char- 
acteristics of Semi- algebraic Sets, Discrete and Computational Geometry, 
22 1-18 (1999). 

[5] S. Basu, Polynomial time algorithm for computing the top Betti numbers 
of Semi-algebraic sets defined by quadratic inequalities, preprint. 

[6] S. Basu, R. Pollack, M.-F. Roy, Computing the Euler-Poincare char- 
acteristic of sign conditions, Computational Complexity, 14:53-71 (2005). 

[7] S. Basu, R. Pollack, M.-F. Roy, Algorithms in Real Algebraic Ge- 
ometry, Springer- Verlag, 2003. 

[8] J. Bochnak, M. Coste, M.-F. Roy, Geometric algebrique reelle, 
Springer- Verlag (1987). Real algebraic geometry, Springer- Verlag (1998). 

[9] G. E. Bredon, Sheaf Theory, Springer- Verlag (1996). 



Computing Euler-Poincare Characteristics 17 



[10] D. Grigor'ev, D.V. Pasechnik, Polynomial time computing over 
quadratic maps I. Sampling in real algberaic sets, Computational Com- 
plexity, 14:20-52 (2005). 

[11] R. M. Hardt, Semi-algebraic Local Triviality in Semi- algebraic Map- 
pings, Am. J. Math. 102, 291-302 (1980). 

[12] J. McCleary A User's Guide to Spectral Sequences, Second Edition 
Cambridge Studies in Advanced Mathematics, 2001. 



Manuscript received March 21, 2005. 

Saugata Basu 
School of Mathematics, 
Georgia Institute of Technology, 
Atlanta, GA 30332, U.S.A. 
saugata@math . gatech . edu 



